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ABSTRACT 



The primary goal of this thesis was the development of an 
experimental technique, and supporting software, for the ac- 
quisition and analysis of unsteady velocity data generated by 
an oscillating airfoil. This research was in support of a 
major investigation of the compressibility effects on dynamic 
sta 1 1 . 

The experimental procedure involved schlieren flow visu- 
alization for comparison of steady and unsteady flow fields, 
and for determination of parameters for further study. Laser 
Doppler velocimetry was employed for obtaining velocity data 
in the airfoil wake. For unsteady data, the airfoil was 
oscillated in pitch about its quarter chord. 

The data analysis produced wake profile plots represent- 
ing the flow field disturbed by the airfoil. Results were 
obtained for steady and unsteady conditions. 
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I . INTRODUCTION 



A. DYNAMIC STALL 

The possible benefits and consequences inherent in 
dynamic stall have stimulated experimental investigations 
leading to significant progress in understanding this complex 
phenomenon. Much of the earlier work in this area was car- 
ried out by McCroskey using a two-dimensional airfoil 
section oscillating in pitch about its quarter-chord . The 
results of this effort succeeded in clarifying the features 
of dynamic stall up to a Mach number of 0.3 CRef. 13. 

Several significant distinctions between static and 
dynamic stall should be observed. Static stall is character- 
ized by massive flow separation and loss of lift coinciding 
with some critical angle of incidence. The aerodynamic 
forces on the airfoil vary uniquely, or nearly so, with angle 
of attack. Dynamic stall, occuring during rapid increases 
in incidence, results in a delay of the associtated loss of 
lift until well beyond the static stall angle. Dynamic stall 
also tends to be considerably more persistant. The flow 
field does not immediately adjust to reductions in incidence 
angle as with static stall. CRef. 23 

Dynamic stall is charac ter ized by the shedding and pas- 
sage of a vortex over the upper surface of an airfoil. The 
production of this vortex begins shortly after the incidence 



Q 



has exceeded that of static stall. 



A thin layer of reversed 



flow in the boundary layer, evolving into a vortex near the 
leading edge, represents the initial phase of a dynamic stall 
cycle. As this disturbance moves across the upper surface of 
the airfoil, large negative pitching moments, significant 
drag rise, and a monotonic increase in lift result. The 
aerodynamic forces, reaching a maximum as the vortex ap- 
proaches the trailing edge, decrease dramatically with the 
shedding of this vortex. Reattachment of the flow will then 
typically occur at some angle less than that corresponding to 
static stall. CRef. 2 U 

B. THESIS GOALS 

The research project generating this thesis represents a 
small part of a major investigation aimed at developing a 
thorough understanding of the effects of compress i b i 1 i ty and 
pitch oscillation on dynamic stall. The primary goal of the 
research associated with this thesis was the development of 
an experimental procedure, and supporting software, for the 
acquisition and analysis of unsteady flows with a laser dop- 
pler velocimetry system. The procedure and software devel- 
oped was to provide an initial analysis of possible areas of 
interest and a base upon which to build as the sophistication 
of the experimental methods increased. A successful attain- 
ment of this goal required obtaining velocity data in the 
wake of an oscillating airfoil and wake profile plots pro- 



duced from this data. 



I I . EXPERIMENTAL CONSIDERATIONS 



A. EXPERIMENTAL METHODS 

1 . Laser DopPler Veloc imetcy: 

Since its introduction in 1964, laser Doppler vel- 
ccimetry has proven to be a highly reliable means for obtain- 
ing measurements in complex turbulent flows. The method is 
capable of performing instantaneous, non- i n tr us i ve velocity 
measurements under a wide range of conditions. Being rela- 
tively independent of fluid properties, accurate measurements 
are possible without the requirement of calibration. CRef. 33 
The basic principle behind laser Doppler velocimetry 
is the measurement of the velocity of particles travelling 
with the fluid. These particles, when illuminated with a 
focused laser beam, become sources of scattered light. The 
Doppler shift of the scattered light is then used to deter- 
mine the velocity of the particles and, depending upon how 
well the particles are following the flow, the fluid velocity 
as we 1 1 , 

As a minimum, six components are necessary for a com- 
plete LDV system^ laser source, transmitting optics, receiv- 
ing optics, photodetector, signal processor, and data 
processor. The function of the transmitting and receiving 
optics are focusing and light collecting, respectively. The 



p hotodetec tor serves to convert light signal to an electrical 



signal which is then used by the signal processor to provide 



a measurement of the frequency. The capability of detecting 
reversed flows requires an additional component for frequency 
shi f ting . 

Due to the wide range of disciplines involved, op- 
tics, electronics, light scattering, and signal processing, 
the implementation of this method can be difficult CRef. 31. 
Consequently, the expertise of the operator can have a pro- 
found impact on the results of LDV measur ements . It is 
therefore imperative, to an even greater degree than many 
other methods, that considerable attention be given to ap- 
plying correct measurement technique. 

2. S cliU^r^n E1S2H V i 1 i S.^t j, gn 

Schlieren flow visualization systems provide a quali- 
tative analysis of the density gradient in compressible flow 
fields. The fundamental concept of this method is the vari- 
ation of the refractive index gradient of a fluid with 
dens i ty . 

A basic schlieren system includes a light source, 
transmitting and receiving mirrors, a knife edge, and a 
focusing screen. The transmitting mirror serves to direct a 
parallel light beam through the test section and onto the 
receiving mirror- An image of the light source is then 
formed at the focal point of this mirror. The knife edge, 
placed in the plane of the light source image, cuts off part 
of the transmitted light and controls the contrast of the 



image focused on the screen. Changes in test section density. 



theref ore , 



produce variations of light intensity and provide 



a visual representat ion of the flow field. CRef. 41 

B. INSTRUMENTATION AND EQUIPMENT 
1 . Equipment Setup 

The experimental procedures for this research were 
conducted at the NASA Ames Fluid Mechanics Laboratory. A 
25-cm by 1 3-cm indraft wind tunnel with a slotted wall test 
section was used. Figure 2.1 shows the test seption with 
plenum chambers above and below the slotted walls. The 
adaptive wall design of this test section was developed for 
use in previous experiments. Although the technique was not 




Figure 2.1 Wind tunnel test section. 



1 3 




used, the slotted walls were helpful in reducing wall 



interference . 

The model for the experiment was a NACA 0012 airfoil 
with 3” chord. The airfoil, attached to the test section 
walls at quarter chord, was connected with mechanical linkage 
to an electric motor, as shown in Figure 2.1. This variable 
speed motor provided a means for oscillating the model about 
a pitch axis and obtaining unsteady data at various 
f requenc i es . 

An optical encoder, also driven by the motor, pro- 
vided 3,600 counts per revolution of 14-bit data on the phase 
angle of the oscillating airfoil. The characteristics of the 
encoder were absolute in that a whole word output, with a 
unique code pattern, represented each of the 3,600 discrete 
positions. This code was derived from independent tracks on 
the encoder disk c or respond i ng to individual photodetectors. 

Test section Mach number and velocity were obtained 
by using a 24-port Scanivalve with a differential pressure 
transducer. Calculation of these values required meas- 
urements from four of the ports, two for Mach number computa- 
tion and two for calibration of the transducer. The linear 
relationship between the measured pressure and the corre- 
sponding voltage produced by the transducer was used for 
determining the pressures . In accomplishing this, one port 
of the Scanivalve, vented to the atmosphere, measured ambient 
pressure and a second port, tapped to a mechanical vacuum. 



pump , 



measured a calibration pressure. 



A conversion factor 



was then established by dividing the difference of the known 
pressures, atmospheric and calibration, by the difference in 
voltages produced by the transducer. 

Location of the LDV system probe volume, or beam 
crossing point, was controlled by a three axis traversing 
mechanism. A Summit Taskmaster numerical control computer, 
originally developed for controlling movements of industrial 
machine tools, was used for this task. The numerical control 
computer NCC ) was capable of manual operation, programmed 
operation, and, with software adaptation, operation from the 
computer workstation. The NCC provided control inputs to two 
sets of electric motors, one for the transmitting optics and 
the other for the receiving optics. The necessity of using 
different traversing motors for the transmitting and receiv- 
ing optics placed a great demand on the synchronization of 
this system. Due to the high sensitivity of the LDV optics 
train, the synchronized movement was often inadequate for 
maintaining signal strength. 

Phase encoder data and digitized velocity data from 
the LDV signal processors were channeled to a 3D-LDV Computer 
Interface (Cl). The Cl, multiplexing the inputs to a single 
digital data channel output, provided 16-bit parallel output 
data to the computer data bus. The Cl was capable of accept- 
ing up to six channels of digital data and from one to three 
event pulses from the signal processors. The multiplexing 
cycle could be initiated by random event pulse inputs or 



synchronized inputs, the latter being necessary for unsteady 



data acquisition. 



When operating in the synchronized mode, 



two or three pulses, as selected, occuring within a set time 
(coincidence time) would be required to initiate the multi- 
plex cycle and feed data to the computer. 

2. Laser Doppler Velocimeter Subsystems 

A two component LDV system was used to obtain the 
flow field velocity data. Figures 2.2 and 2.3 show the gen- 
eral arrangement of the system components for the transmit- 

I 

ting and receiving optics sides, respectively. 

A 5 Watt Argon ion laser, Spectra Physics Model 164, 
was used for the system. The unit was operated in the 




Figure 2.2 LDV transmitting optics 




Figure 2.3 LDV receiving optics. 

multiline mode and produced a vertically polarized beam com- 
posed of all wavelengths inherent to this type of laser. 

Figure 2.4 illustrates the sequence of steps required 
of the transmitting optics. A set of mirrors was first used 
to lower the multiline beam to the level of the optical 
components. After passage through a' collimatorf used to 
ensure collimation of the beam and subsequent crossing of the 
beams at their waistSf the beauns pass through an attenuator. 
The function of this component was to reduce the beam intens- 
ity without inhibiting transmission of any wavelengthsi a 
possible consequence of simply reducing the power output of 
the laser. Prior to separation of the wavelengths by the 
dispersion prismi it was necessary to horizontally polarize 
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the bean by passing it through a half-wave plate, 



or po 1 ar i - 



zation rotator, in order to obtain maximun light transmis- 
sion. The two stronger wavelengths, 514.5 nm (green) and 486 
nm (blue), were then guided by turning mirrors to the remain- 
ing optical components. 

Before reaching their respective Bragg Cells, the 
green beam was split vertically and the blue was displaced 
to the center of the component assembly, vertically polar- 
ized, and split horizontally. The relationship between the 
splitting and polarizing planes was driven by signal effi- 
ciency. High si gnal-to-noise ratio requires split beams to 
have equal intensities. The beams will be of approximately 
equal intensity when the polarity of the incident beam is 
perpendicular to the plane of the exiting beams CRef. 3]. 
Therefore, the polarity of the vertically split beam was 
maintained as horizontal and that of the horizontally split 
beam was rotated to vertical. As later discussed, this dif- 
ference in polarity was also necessary for separation of the 
two wavelengths by the polarization beamsplitter in the re- 
ceiving optics. 

Bragg Cells are acousto-optic devices used to gen- 
erate a frequency shift in the beam and allow for resolution 
of the 180 degree directional ambiguity. One beam of each of 
the two wavelengths was directed through a Bragg Cell pro- 
ducing a 40 MHz shift. The unshifted beams were passed 
through glass rods to compensate for the path length differ- 
ence experienced by the shifted beam. Due to the nature of 



the Bragg Cell, the f requenc y- sh i f ted beam was deflected from 
its initial alignment. Therefore, beam alignment prisms, or 
steering wedges, within the Bragg Cell modules served to 
bring the beams back to parallel. 

A microscope objective, placed in the test section at 
the probe volume, was used to observe the beams as adjust- 
ments were made to ensure the beams focused to a common 
point. If the two probe volumes are not at the same location 
within the test section, then the measurements will be made 
at different points. The necessity for all four beams to 
overlap at the focal point was especially prevalent for the 
acquisition of unsteady data due to the need for a high 
coincidence data rate. 

Receiving optics are designed to image, to the great- 
est extent possible, the light scattered at the probe volume 
onto the photodetectors. Since the scattered light intensity 
is the greatest in the forward direction, a forward scatter 
system is most often preferred for obtaining the greatest 
signal-to-noise ratio. The system used for this analysis was 
an off-axis, forward scatter configuration. Off-axis col- 
lection angles usually allow good noise isolation by admit- 
ting only light scattered at the sample volume and improve 
measurement resolution by reducing the probe volume length. 
[Ref. ?1 

The first three lenses in the receiving optics were 
used to collect scattered light from the probe volume and 
collimate the beam (Figure 2.5). The beam was then directed 



20 



to a polarization beamsplitter for separation of the two 
colors. The two components were then passed through line 
filters, for blocking ambient light and preventing cross 
talk, and to lenses used for focusing the beams onto the pin- 
hole aperature of their respective photomultiplier tubes. 
The photomultiplier tubes converted the light energy into 
electronic signals which were then passed to signal proces- 
sing equipment. Pr eamp 1 i f i ers were located in close proxim- 
ity to the photomultiplier tubes and used to increase the 
signal level bef ore the i n tr oduc tion of additional noise into 
the cables. 

Signal processing was accomplished with Macrodyne 
counter processors. Counter processors operate by employing 
a crystal controlled oscillator to time over a selected num- 
ber of cycles within the burst signal. The clock count is 
initiated when the filtered signal exceeds a set threshold. 
The Doppler frequency is then determined by the time required 
to count a specified number of cycles. The Macrodyne proces- 
sors use a 125 MHz clock and count over 8 or 16 cycles, as 
se 1 ec ted . 

The Macrodyne processors use a periodicity check 
coupled with a three-level validation scheme to reject noise. 
The periodicity check uses two counters to make period meas- 
urements over a different number of cycles, 5 and 8 in this 
procedure. If the ratio of the two measurements is not with- 
in a specified percentage of 5/8, the signal is rejected. 
Three-level validation requires a signal to exceed a positive 
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threshold level 



pass through zero, and then exceed a nega- 



tive threshold level before the next zero crossing. If these 
conditions are not met on any cycle the counter will reset 
and immediately process a new Doppler burst. This method 
prevents rejection of potentially good signals having noisy 
cycles at the beginning of the burst. 

3 . Schl ieren System 

The Schlieren system included two concave mirrors, 
each with a 120 inch focal length, a stroboscopic light 
source, a vertical knife edge, a focusing lens, and a Polar- 
oid film support for obtaining photographs of the flow field. 
The stroboscope, a type 3015 Strobrite, has a maximum output 
intensity of 300 million beam candles and, using an internal 
oscillator, can provide flashing rates in excess of 1,000 per 
second. By adapting a lens cover with a vertical slit to the 
stroboscope a slit source was obtained and the resolution of 
the image significantly improved. 

To obtain the unsteady photographs a digital compar- 
ator was used to trigger the stroboscope and provide the 
light for the photographic imaging on the Polaroid film. The 
comparator received 14-bit phase input from the optical en- 
coder and provided the triggering pulse when the 
instantaneous phase angle of the airfoil coincided with the 
angle specified to the comparator. 
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DATA ACQUISITION AUTOMATION 



C . 



A Digital VAXstation II minic om puter with graphics work- 
station was used to automate the data acquisition process. 
The signal path is illustrated in Figure 2.5 starting from 
the probe volume, through the signal processing phases, and 
terminating within computer memory. The main program used in 
accomplishing this task, *LDV_ACQ*, was originally developed 
to read data f rom the analog to digital converter. This was 



later modified to obtain pressure data for calculation of the 



test section velocity and call supporting subroutines for the 
acquisition and analysis of the velocity data. 

The subroutines assisting with various aspects of the 
acquisition process are listed below and included in the 
appendix of this report. 

(1 ) HWTRVT - allows control of the traversing mechmisn 
from the computer terminal by sending character data 
to the interface. 

(2) HWWA - combined with an internally called subrou- 
tine, DAMUX , retrieves data from multiplexer, inter- 
prets phase and velocity information, creates data 
files, and generates histogram arrays and scaling 
parameters . 



( 3 ) 



PT_PLOT - 
veloc i ty 
angles . 



provides graphical display of computed 
values correspond i ng to airfoil phase 



(4) PPLOT - similar to PT_PLOT but provides the addit- 
ional capability of allowing the user to vary the 
plotted velocity range. 



( 5 ) 



HIST - plots 
points which 



a histogram of the percentage 
fall within specified velocity 



of data 
limits. 
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III. QAJh ACQUISITION 



A. TEST SECTION VELOCITY 

The main program for the acquisition process, *LDV_ACQ’, 
was responsible for the calculation of test section velocity. 
In combination with scanivalve positioning subroutines, the 
program acquired pressure transducer data from 7 ports. As 
previously discussed, measurements from 2 of the ports were 
used for calibration. Pressure data from 3 ports was dis- 
carded, and the other two provided the static pressure at the 
entrance to the test section and the total pressure measured 
at the nozzle inlet. Assuming isentropic, perfect gas condi- 
tions, this supplied the necessary parameters for the Mach 
number computation by employing the following relationship: 

with a value of 1.4 used for 7 • T'he test section velocity 
was then calculated by: 

V = M X a 

where a, the sonic velocity, was obtained as follows: 

a = \/ yRT 

The temperature required for this calculation was read from a 
wall thermometer in the test cell and entered when prompted 
by the program. 

The velocity values were transferred to data files and 
later used to normalize the LDV velocity data. 
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B. SCHLIEREN PHOTOGRAPHY 



A qualitative analysis of the steady and unsteady flow 
fields was conducted at two Mach numbers, 0.3 and 0.5, and 
angles of incidence varying from 0 to 20 degrees. In ac- 
complishing this, test section Mach number was set by adjust- 
ing the wind tunnel throat and verified with Scanivalve 
pressures and calculations f rom *LDV_ACQ’. Final adjustments 
to the Stoboscope intensity and knife edge were then made to 
ensure optimum contrast. 

C. WAKE VELOCITIES 

With the desired wind tunnel conditions established and 
the probe volume in position for the data point, ’HWWA’, a 
subroutine called by the main program , was employed for raw 
data acquisition and conversion to velocity and phase infor- 
mation. An internally called subroutine, 'DAMUX*, acquired 4 
or 6 channels of data, as selected, from the LDV Computer 
Interface (Cl). The Cl was capable of transferring either 2, 
4, 6, or 8 channels of data, with the first two of the 
channels reserved for time and status information generated 
by the Cl. During the acquisition of steady state velocities, 
4 channels were selected to provide the 2 velocity components 
in addition to time and status data, which were discarded 
since they were not relevant for the experiment. For unsteady 
acquisition, 6 channels of data were transferred to accom- 
modate the additional data channel required for phase input. 
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As previously mentioned, the signal processors determined 
the Doppler frequency by using the time required to count a 
specified number of cyles, 8 for this analysis. This raw 
data was transferred to the computer in 16 bits, with a 10 
bit mantissa, a 4 bit exponent, and 2 nonessential bits used 
for controlling data flow. Due to the format of the time 
words, it was necessary to invert and reverse the bits before 
extracting the exponent and mantissa values. The signal 
period was then calculated as follows: 

ts = ( 1 /32 )( D9D8. . .DO ) x 2 

where D9D8...D0 represented the binary mantissa and EXP was 
the range selected with the processor, used to divide the 
number of clock counts in the Doppler burst. 

The signal frequency, the inverse of the period, was ob- 
tained from: 

= (32 X 1000 )/(2^^^x mantissa) 

The signal frequency was composed of the Doppler frequen- 
cy, the 40 MHz shift frequency, and the downmixing frequency. 
This allowed the Doppler frequency to be deduced as follows: 

f =f -f , 

4 I ahlff mil 

The 8 cycles counted within the Doppler signal physically 
represent 8 fringes crossed by a particle travelling with the 
flow. Consequently, the velocity component perpendicular to 
the fringe pattern may be calculated by the product of the 
Doppler frequency and the fringe spacing. The equation for 
the latter is: 
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<5 = X / [2 • sin(d/2)] 

where X is the wavelength of the beam and 9 is the angle be- 
tween the intersecting beams. 

After obtaining a sufficient number of data points at a 
probe position (2,000 for unsteady flows), the velocity and 
phase data, as well as the probe position, test section 
velocity, and mixing frequencies, were stored in files for 
later analysis. The Taskmaster traverse controller was then 

I 

used to move the probe volume and data collection continued- 
The data acquisition software included two routines for 
on-line performance checks of the data. ’HIST* was used to 
provide a histogram of the percentage of data points which 
occurred within specified velocity limits. Figure 3.1 shows 
a sample plot provided by the subroutine. This data valida- 
tion method is most useful for steady flow conditions, in 
which case a high percentage of data concentrated around the 
estimated velocity would be a positive indicator. For un- 
steady analysis ’PT_PL0T* provided a graph of computed ve- 
locity values corresponding to airfoil phase angles. With 
this information available, the data quality could be 
easily assessed and a determination made regarding the ac- 
ceptance of these values. 
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PERCENTAGE OF DATA 



VERTICAL VELOCITY HISTOGRAM 




VELOCITIES 

a 



Figure 3.1 Sample histogram. 
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IV. BSSV LTg AJIC DISCUSSION 



A. SCHLIEREN PHOTOGRAPHS 



The Schlieren photographs obtained from the analysis of 
the steady flow field were intended to provide information on 
the parameters corresponding to early flow separation and a 
reference for comparison with unsteady flows at, the same 
conditions. Steady state separation was first observed at a 
Mach number of 0-5 and a 10 degree angle of incidence (Fig- 
ure 4.1 ). Based on this, the parameters for the unsteady 




Figure 4.1 Steady state Schlieren photograph; 
mean a =1 0 degrees. 
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analysis were selected to be a Mach number of 0.5 and a phase 
angle corresponding to a 1 0 degree angle of incidence. Un- 
steady flow photographs were then obtained at oscillation 
frequencies of 5 and 25 Hz for comparison with each other and 
with wake profiles generated from LDV data (Figures 4.2 and 
4.3). 

From a comparison of the steady and unsteady photographs, 
it would appear that oscillating the airfoil reduces the size 
of the separated wake and, consequently, mitigates the effect 
of stalling. This is more apparent at the higher frequency 




Figure 4.2 Unsteady flow Schlieren photograph; 
f=5 Hz, a =10 degrees. 
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Figure 4.3 Unsteady flow Schlieren photograph; 
f=25 Hz, a =1 0 degrees. 



oscillation, a logical result due to the greater effect of 
time unsteadiness at higher reduced frequencies. 

B. DATA ANALYSIS SOFTWARE 

The main program, •ORGIZE*, and its supporting subrou- 
tines were designed to provide wake profile plots from the 
data files generated during acquisition. To accomplish this, 
the data file corresponding to each point in the profile is 
input to the software. A series of routines first reorgan- 
izes the data files and discards statistically bad data 
points, those outside of two standard deviations from the 
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mean. Due to the oscillatory nature of the forcing function 
producing the unsteady data, a sinusoidal curve fit is then 
introduced to allow evaluation of the first harmonic of the 
velocity distribution over the phase. A unique sinusoid is 
generated for each data point within the profile. 

A minor degree of uncertainty is introduced into the ve- 
locities by the data interpolation and curve fitting rou- 
tines discussed above. However, unlike many other velocity 
measurement techniques, LDV data is obtained randomly due to 
the random arrival of the seed particles at the prove volume. 
Without this curve fit, it would be difficult to generate 
velocity information at the same desired phase angle for each 
data point in the profile. 

During the generation of the wake profile the user is 
provided a comparison plot of mean values and the sinusoid 
curve fit. The mean values are provided by ’BDPTS’ and com- 
puted within each 20 degree phase window, as illustrated in 
Figure 4.4. This was helpful in ascertaining the quality of 
the curve fit and, consequently, the quality of the wake pro- 
file results . 

The data analysis routines are summarized below and in- 
cluded in the appendix of this report. 

(1) ORGIZE - main program , calls supporting subroutines; 
sequentially organizes data to assist with subse- 
quent analysis. 

(2) PT_PL0T - provides graph of computed velocity values 
correspond! ng to phase. 

(3) BDPTS - computes mean and standard deviation of ve- 
locities occurring within 20 degree phase windows; 
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VELOCITY (ft/s) 



HORIZONTAL VELOCITIES 




discards data points outside of two standard devia- 
tions of the mean within each of these windows; pro- 
vides standard deviation plots of data. 

(4) AVC - averages data valued occurring at same dis- 
crete phase. 

(5) INTER? - provides a linear interpolation between 
data, resulting in 180 equally spaced points; neces- 
sary for input to FOURIER. 

(6) FOURIER - computes coefficients for a sinusoidal 
curve fit of the data adapted from June 1 977 OCS 
Airfoil program ) . 

(7) CVFIT - fits sinusoidal curve to velocity data ( u- 
tilizes FOURIER). 
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C8) PLOT - provides graphical comparison of curve fit 
with mean values generated by BDPTS. 

(9) PROPLOT - provides wake profile plots. 

C. WAKE PROFILES 

LDV data was obtained for steady and unsteady flows at a 
Mach number of 0.5 and a 10 degree angle of incidence. The 
steady profiles were collected at distances of 0.5, 1.0, and 
2.0 inches behind the trailing edge of the airfoil. These 
values were non-d imens iona 1 i zed by the chord length of the 
airfoil and are indicated as such on the profile plots. Due 
to time limitations, only one data profile was collected for 
the unsteady case, at two inches behind the trailing edge. 

The wake profile plots are presented in the following 
pages of this section. A lack of data points in the lower 
part of the flow field will be apparent upon inspection of 
these results. This was due to a mechanical restriction on 
the traversing mechanism. 

The curve fit of the plotted values of the steady data 
was approximated and hand drawn to aid the reader in ob- 
serving the effects of the airfoil in the flow field. The 
software provides a spline interpolation fit between the 
data, as used in the unsteady plots. However, due to the 
amount of data scatter involved, this method was not used for 
these results. 



35 



HORIZONTAL VELOCITY PROFILE 




• VELOCITY RATIO (U/UdNFIN)) 



Figure 4.5 Steady wake profile. 
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PROBE POSITION (IN. FROI CENTER! 



VERTICAL VELOCITY PROFILE 




VELOCITY RATIO (V/U(INFIN)) 



Figure 4.6 Steady wake profile. 
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PROBE POSITION (IN. FROM CENTER) 



HORIZONTAL VELOCITY PROFILE 




VELOCITY RATIO (U/U(INFIN)) 



Figure 4.7 Steady wake profile. 
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PROBE POSITION (IN. FRijM CENTER) 



VERTICAL VELOCITY PROFILE 




VELIXITY RATIOl(V/U(ITfIN)) 



Figure 4.8 Steady wake profile. 
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PROBE POSITION (IN. FROM CENIERI 



HORIZONTAL VELOCITY PROFILE 




VELOCITY RATIO -(U/UdNFIN)) 



Figure 4.9 



Steady wake profile. 
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PROBE POSITION (IN. FROM CENTER) 



VERTICAL VELOCITY PROFILE 




VELOCITY RATIO; (V/U(INFIN)) 



Figure 4.10 Steady wake profile. 
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PROBE POSmON (IN. FRlJM CENIER) 



HORIZONTAL VELOCITY PROFILE 
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Figure 4.11 Unsteady wake profile. 
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PROBE POSITION (IN. FROM CENIER) 



VERTICAL VELOCITY PROFILE 




VELOCITY RATIO. (V/U(INFIN)) 



Figure 4.12 Unsteady wake profile. 
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V. CONCLUSIONS AND RECOMMENDATIONS 



The purpose of the results obtained during the experi- 
mental phase of this thesis was to validate the technique and 
software developed. This was accomplished when unsteady 
data, generated by an oscillating airfoil, was successfully 
acquired, analyzed, and reduced to wake profile plots. 

Throughout the course of this investigation hardware 
problems have a posed a signigicant challenge to acquiring 
data within the existing time constraints. Consequently, the 
data obtained was accepted and included in this report with- 
out the opportunity to more carefully investigate regions of 
significant data scatter and improve, or validate, the data 
plots . 

The results produced by these procedures confirm a sound 
technique, but also reveal several areas in need of addi- 
tional emphasis. At present, an interpolation subroutine is 
required to produce data at equally spaced intervals before 
generation of the curve fitting coefficients by ’FOURIER’. 
Combined with the inability to utilize a large number of data 
points ( < 300), this limitation will inherently introduce 
additional uncertainty into the results. A second software 
modification to be considered is the generation of the veloc- 
ity data at various phase angles simultaneously. Currently, 
the data must be processed for each separate phase angle 
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desired. This would eliminate the need for storage of large 
amounts of raw data and significantly reduce the processing 
time. By increasing the range of vertical movement of the 
traverse, the mechanical restriction that prevented the ac- 
quisition of data in the lower regions of the wake can be 
el iminated . 

Continuation of this study will be accomplished with a 
newly developed test section currently being installed. This 
test section will enable schlieren flow visualization and LDV 
data acquisition at locations over the airfoil which were 
previously blocked by the mechanical linkage used for sup- 
porting and oscillating the airfoil. This, combined with the 
modifications suggested above, will enable an in-depth analy- 
sis of wake and boundary layer characteristics associated 
with the dynamic stall phenomenon. 
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APPENDIX 



COMPUTER PROGRAMS 



PROGRAM LDV ACQ 



Adapted froa LABSTAR Exasple3_l by S. D. Hedrick in Hay, 1987 for reading 
A-D data. Modified in August? 1987 to be used as the central program in 
the acquisition of unsteady velocity coaponents with an LDV. 



INCLUDE 'sys$library:LIOSET.rOR' 

C declare local variables 

INTEGER sys_stat Istatus ret by LIO calls 

INTEGER devoid, clk_id ILIO device ID 
INTEGER data_length~ tnuaber of bytes of data read 
INTEGER i, npt, dnpt 
INTEGER chn, nsteps, ier 

REAL rate, volt(0:23), press(0:23), corfac 

REAL SUM, aean, calp, conv, aab, teap, aach, velo 

CHARACTER*! ans 

CHARACTER*4 clock 

LOGICAL hoae 



C declare data buffers 

INTEGER*2 raw_data( 501 ) 
REAL VOLTAGES ( 501 ) 



OPEN(UNIT-13, FILE-' scani .dat ' , STATUS- ' NEW' ) 

write(*,*) 'ENTER A-D CHANNEL TO BE USED (0,1, 8,9)' 
read! * , * ) chn 

10 write(*,*) 'ENTER CLOCK RATE ( . 01 - 4500 Hz)— ALLOW FOR A', 

1 ' MAXIMUM OF 2 SECONDS TO ACQUIRE DESIRED SAMPLES' 

read( * , * ) rate 

if (rate .It. .01 .or. rate .gt. 4500) then 

write(*,*) 'RATE IS NOT WITHIN SPECIFIED LIMITS' 
go to 10 
end i f 

20 write(*,*) 'ENTER NUMBER OF DATA POINTS TO BE COLLECTED', 

1 ' ( 500 aax) ' 

read(*,*) npt 

if (npt .It. 1 .or. npt .gt. 500) then 

write(*,*) 'NUMBER OR POINTS IS NOT WITHIN LIMITS' 
go to 20 
end if 

npt - npt ♦ 1 
dnpt - npt * 2 

write(*,*) 'ENTER CLOCK DEVICE (K2AO,K2BO)' 
read( * , ' ( A4 ) ' ) clock 

Attach the ADVll-D and set up for queued I/O 



sys_Stat-LIO$ATTACH(dev_id, 'A2A0' , LIO$K QIO) 

IF(.NOT. sys_stat) CALL LIB$SIGNAL( %val ( sys_sta t ) ) 



Set up the ADV for synchronous transfer on selected channel 
Gain of 1 



sys_stat - LIO$SET_I (dev id, LIO$K SYNCH, 0) 

IF(.NOT. sys_statT CALL LIb5siGNAL( %val ( sys_stat ) ) 
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sys_stat - LIO$SET_I (dev id, LIO$K AD^CHAN, 1, chn) 

IF(.NOT. sys_statT CALL LIB^S IGNAL( %val ( sys^stat ) ) 

£ys_stat - LIO$SCT_I (dev id, LIO$K AD^GAIN, 1, 1) 

IF( .NOT. sys^StatT CALL LIB^S IGNAL( Ival ( sys_stat ) ) 

Get a raw_data buffer 

length of buffer and returned data length is in bytes 
Attach clock and set rate 

ho*e - .true, 
nsteps - 0 
do 40 j *0 , 7 

call dalljllab(home, nsteps, ier) 

sys stat- LIO$ATTACH( clk^id, clock, LIO$K_QIO) 

“ IF(.NOT. sys stat) CALL LIB$S IGNAL( Ival ( sys^stat ) ) 
sys stat-LIO$SCT_R( cTk_id, LIO$K_CLK RATE, l,rate) “ 

IF(.NOT. sys stat) CALL LIB5 s IGNAL( % val ( sys_stat ) ) 
sys Stat-LIO$SET_I( cIk_id,LIO$K_TRIG, 1, LIO$K IMMEDIATE) 

IF( .NOT. sys stat) CALL LIB$SIGNAL( %vaT( sys^stat ) ) 
sys_stat • LIO$READ(3ev id, raw^data , dnpt ,data length, ) 

IF( .NOT. sys stat) CALL LIB$S IGNAL( %vaT( sys Stat)) 
sys_stat - LIO$SET l“(dev_id, LIO$K_TRIG, 2, LIO$R_^LK_POINT 
1 clk_idT 

IF( .NOT. sys_stat) CALL LIB$S IGNAL ( I val ( sys_s tat ) ) 
Convert the raw data to voltages 

call LSP$FORMAT_TRANSLATE_ADC( raw^dat a, voltage s,npt, , , ) 



calculate the mean value of the voltages and detach the clock 

sum • 0.0 
do 30 i-2,npt 

sum • sum voltages(i) 

30 continue 

mean - sum / float(npt-l) 
vol t ( j ) - mean 

sys_stat - LIO$DETACH( elk id, ) 

IF( .NOT. sys_statT CALL LIB$SIGNAL( %val ( sys_stat ) ) 

40 continue 

detach from the A/D 

sys stat - LIO$DETACH(dev id, ) 

IF(.NOT. sys_statT CALL LIB$SIGNAL( %val ( sys_sta t ) ) 

convert scanivalve voltages into units of pressure 

write(*,*) 'ENTER CALIBRATION PRESSURE (in. of Hg ) ' 
read( * , * ) calp 

write(*,*) 'ENTER REQUIRED CORRECTION FACTOR FOR CALIBRATION 
1 ' PRESSURE ( in. of Hg) ' 

read(*,*) corfac 
calp-calp-fcorfac 

wrltt(*,*) 'ENTER ATMOSPHERIC PRESSURE (in. of Hg)' 
r«ad(*,*) amb 
praas(O) - calp 
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press(l) • amb 

conv - ( pres6( 1 )-press ( 0 ) ) / ( vol t ( 1 ) -vol t ( 0 ) ) 
do 50 j-2,7 

press(j) • conv * volt(j) calp - volt(O) * conv 

50 continue 

WRITE( * , * ) VOLT( 0 ) , VOLK 1 ) , VOLT( 2 ) ,VOLT( 7 ) 
write(13,55) 

55 format (2X, 'PORT' ,4X, 'PRESSURE' ,//) 

do 60 i *0 , 7 

write( 13, ' ( 3X, 12, 5X, F9. 4) ' ) i, prc6S(i) 

60 continue 

compute test section velocity 

write( *,*) 'ENTER ANBIENT TEMPERATURE (deg. F)' 
read(*,*) temp 

mach-sqr t ( abs (5.*( (pres6(2)/pres6(7) )**(. 4/1 . 4 )-l . ) ) ) 
velo-aach* sq r t (1.4*53. 3*32. 174*( temp'*>4 59 . 67 ) ) 

write(*,*) 'TEST SECTION MACH NUMBER AND VELOCITY — ',mach,velo 

write(*,*) 'CONTINUE WITH LDV DATA ACQUISITION? (Y/N, DEF-Y)' 
read( * , ' ( A1 ) ' ) ans 

if (ans ,eq. 'n' .or. ans .eq. 'N' ) go to 70 

call subroutine to move LDV traversing mechanism and acquire data 

call HWTRVT(velo) 

C 

70 stop 
end 
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SUBROUTINE HWWA ( XPOS , Z POS , VELO ) 

ACRONYM: Hardware checkout - test of DRVll-WA 

PURPOSE: 

For LVIS, HWDVYB drives the acquisition of the LDV data. 

C It is a test program for the YB driver. 

C 

C METHOD: 

C 

C BEGIN 

C Install loadable drivers if necessary. 

C Read the various instrument types in succession, skipping those which 
C are not to be read. 

C Remove loadable drivers if necessary. 

C END 
C 

C PARAMETERS: 

C ARC DIM TYPE I/O/S DESCRIPTION 

C 

C COMMONS USED: 

C 

C FILES USED: None 

C 

C ERROR HANDLING: 

C lER - 1 means no "fatal" error. 

C 

C NOTES: 

C 

C LOCAL VARIABLES; 

C VAR DIM TYPE DESCRIPTION 

C 

C EXTERNAL REFERENCES: 

C NAME DESCRIPTION AND SOURCE 

C DVWA Calls DRVll-WA device driver to acquire A/D data. 

C 

C STANDARDS VIOLATIONS: None. 

C 

C ENVIRONMENT: DEC VAX/VMS AND FORTRAN 77 

C 

C DEVELOPMENT HISTORY: 

C DATE INITIALS 

C 07/26/83 TML 

C 12/12/85 CLH 

C 

C 02/25/87 GBG 

C 05/06/87 CLH 

C 06/08/87 RRR 

C JUL-AUG,1987 SDH 
C 
C 
C 
C 
C 
C 

C AUTHOR(S): Ted Lichtei 
C 

C. ********* *******«#••* 

c 

IMPLICIT NONE 



DESCRIPTION 
Adapted from DAPNT 

Modified to print results in either 
octal or integer. 

Modified for use with DRVll-WA and VMS 
Modified to accept "new" DAMUXl . 

Modified to subroutine ♦ deleted CLH 5/87 Mod 
Modified to calculate velocities using variable 
mixing frequency as user input, interpret phase 
angle information generated by digital encoder, 
transfer test section velocity and probe posit- 
ion to data file, and generate histogram arrays 
and scaling parameters 



nstein. Informatics General Corp. 
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c 



C 



C 

c 

100 



PARAMETER 

PARAMETER 

PARAMETER 



NRAWMX-80000 

LUNTI-5 

LUNOUT-6 



INTEGER*2 
INTEGER*2 
INTEGER*2 
INTEGER 
INTEGER 
INTEGER 
INTEGER 
INTEGER 
INTEGER 
REAL 
LOGICAL 
CHARACTER*10 
CHARACTER*! 
REAL *4 
REAL*4 
INTEGER 
REAL 



IRAWDT(NRAWMX) , IfORM, MUXTIME (16384) 

MUXSTS ( 16384 ) ,RAWDA( 20000 ) ,RAWDB( 20000 ) 
RAWDC( 20000) 

STATUS, I ER, PHAN 

NCOLUM, lYESNO, NRAWIN,I,J 

LIBSSUBX, NCHANS, NSAMPLES, INDEX, JNDEX 

MRAWA , MRAWB , MRAWC , MRAWD 

I RAWDTA , I RAWDTB , I RAWDTC , I RAWDTD 

IRAWDTE, IRAWDTF 

FRE0(2) ,VEL(2) ,MIXA,MIXB,XP0S,2P0S, VELO 
DATSELCT (3) 

INFILE 

ANS 

YU( 160 ) , YV( 120) , YUK 160 ) , YVl ( 1 20 ) , YMAXU , YMAXV 
SAVE 

COUNT, XLOW,XHI 
X( 4000) ,Y( 4000) 



INCLUDE ' (SSYSSRVNAM) ' 



WRITE(*,*) 'ENTER DESIRED NAME Of DATA FILE— E.G., LDVXXX . DAT , ' , 

1 ' WHERE THE "XXX* MAY BE USED AS A SERIAL ID FOR THE', 

1 ' PROBE POSITION' 

READ( * , ' ( AlO ) ' ) INFILE 

open ( uni t-11 , fi le«INFILE, status* ' new' ) 



CONTINUE 

C * Fill the raw data buffer with this bit pattern: 
C * "1010101010101010" 

C * This equals 125252(8), -21846(10). 

C 

DO 200 I - 1, NRAWMX 

IRAWDT(I) - '125252'0 
200 CONTINUE 



DO I - 1, 16384 

MUXTIME (I) - '125252'0 
MUXSTS (I) - '125252'0 



END DO 



C 

220 CONTINUE 

* Proapt for no. of words to acquire: 



1 

1 



C 



WRITE(*,*)' ' 

WRITE( *,*) 'Starting LDV acquisition' 

WRITE(*,*)' ' 

WRITE( *,*) 'ENTER NUMBER OF CHANNELS TO BE USED', 

' (assuaes non-oscillating airfoil if', 
' less than six)' 

REAO(*,*) NCHAN5 
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224 

225 



226 



C 

230 



9 



10 

C 



WRITE (•, 225) 

rORHAT (' Enter number of samples:') 

READ (*, •) NSAWPLES 

ir (NSAHPLES .GT. 20000) THEN 
WRITE (*, 226) 

FORKAT (' The number of samples must be 20,000 or less.') 
GO TO 224 

END IF 



CONTINUE 

Prompt for output format of data: 

I FORM- 2 

IF ( IFORM.lt. 1 .OR. I FORM. GT. 2 ) GO TO 2 30 



Select the data to be returned. 



DATSELCT ( 1 ) - . false . 

DATSELCT (2) - .false. 

DATSELCT (3) - .TRUE. 

Prompt for the mixing frequencies selected by user 

WRITE( *,♦)' ENTER MIXING FREQUENCIES FOR EACH COMPONENT - (blue, 
Igreen) - in MHz' 

READ(^,*) MIXA,MIXB 



Read and display the data one sample at a time. 

IF (NCHANS .LT. 6) THEN 

WRITE(*,*) 'ENTER AIRFOIL PHASE ANGLE (tenths of degrees, 14)' 
READ(*,*) PHAN 
DO 9 I-1,NSAMPLES 

CALL DAMUX ( NCHANS , 1 , DATSELCT , MUXTI ME , MUXSTS , I RAWDT , I ER ) 

RAWDA( I )-IRAWDT( 1 ) 

RAWDB( I )-IRAWDT( 2 ) 

RAWDC( I )-PHAN 
CONTINUE 
GO TO 15 
END IF 

DO 10 I-1,NSAMPLES 

CALL DAMUX (NCHANS, 1, DATSELCT, MUXTIME, MUXSTS, 
I RAWDT, lER) 

RAWDA( I )-IRAWDT( 1 ) 

RAWDB( I )-IRAWDT( 2 ) 

RAWDC(I)-IRAWDT(3) 

CONTINUE 

DO 19 t-l,NSAMPLES 
HRAHA-O 
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HRAWB-0 

NRAWC-0 

NRAWO-0 

MRAWA-IIBITS( RAWDC( I ) , 12 , 2 ) 

NRAHB-IIBITS( RAWDC( I ) ,8, 4 ) 

NRAHC-IIBITS( RAWDC( I ) ,4 ,4 ) 

MRAWD-IIBITS( RAWOC( I ) , 0, 4 ) 

RAHDC( I )-HRAWA*1000+NRAWB*100+HRAWC*10+HRAWD 

19 CONTINUE 

15 DO 20 I-1,NSAHPLCS 

RAWDA( I )-INOT( RAWDA( I ) ) 

RAHDB( I )-INOT( RAHDB( 1 ) ) 

IRAWDTA-0 

IRAWDTB-0 

IRAWDTC-0 

IRAWDTD-0 

IRAHDTE-0 

IRAWDXr-0 

IRAWDTA-IIBITS( RAWDA( I ) , 10, 4 ) 

IRAHDTB-IIBITS( RAWOA( I ) ,0,10) 

I RAWOTC- 2 * * I RAWDTA • I RAWDTB 
ir (IRAWDTC .EQ. 0) THEN 
WRITE( •, ♦ ) ' IRAWDTC-0' 

GO TO 20 
END IF 

FREOd )-( (32.*10.»*3.)/( FLOAT( IRAWDTC) ) )-40.+HIXA 
I RAWDTD- 1 1 B I TS ( RAWDB ( I ) , 1 0 , 4 ) 

I RAWDTE- 1 1 B I TS ( RAWDB ( I ) , 0 , 1 0 ) 

' IRAWDTF-2**IRAWDTD»IRAWDTE 
IF (IRAWDTF .EQ. 0) THEN 
WRITE( • , * ) ' IRAWDTF-0' 

^ GO TO 20 
END IF 

FREQ(2 )-( ( 32 . *10.**3. )/( FLOAT ( IRAWDTF) ) )-40.+MIXB 
vel( 1 )-Ireq<2 )*0.5145/( 2.*sin( atan( 13./(2.*482.6 ) ) ) ) 
vel(2)-£req( 1 ) *0 . 4880/( 2 . »sin( atan( 1 3 ./( 2 . *482 . 6 ) ) ) ) 
veia)-vel(l)/.3048 
vel(2 )-vel( 2 )/. 3048 

WRITE! 11 ,700) vel(2) , vel( 1 ) ,RAWDC( I ) 

20 CONTINUE 
WRITE! 11 , 900) XPOS,ZPOS 
WRITE! 11 , 950) VELO.NSAHPLES 
WRITE(11,975) HIXA,HIXB 



WRITE! *,*) 'DO YOU DESIRE A HISTOGRAM OF THE VELOCITIES OBTAINED?', 
1 ' I Y/N - DEF-N) ' 

READ! • , ' ( A1 ) ' ) ANS 

IF (ANS .NE. 'y' .AND. ANS .NE. 'Y') THEN 
GO TO 35 
END IF 

CLOSE!UNIT-ll ) 

OPENlUNIT-ll .FILE-INFILE, STATUS-'OLD' ) 

DO 40 I-1,NSAMPLES 

RBAD(11,800) X(I),Y(I) 

40 CONTINUE 
XLOW--200 
XHI--195 
DO 60 J-1,160 



52 



o o o o o o o 



COUNT-0 

DO 50 I-l,NSArtPLES 

IF (X(I) .GE. XLOW .AND. X(I) .LT. XHI ) THEN 
COUNT-COUNT+1 
END IF 
50 CONTINUE 

YU( J)-( FLOAT(COUNT)/NSANPLES)*100 . 

XLOW-XLOW^5 
XHI-XHI^5 
60 CONTINUE 

Calculating a scaling parameter for the horizontal velocity histogram 

DO 62 1-1,160 
YUK I )-YU( I ) 

62 CONTINUE 
DO 63 1-1,160 

DO 65 J-1,159 

IF (YUKJ^l) .LT. YUl(J)) THEN 
SAVE-YUK J) 

YUK J)-YU1( J^l ) 

YUK J-fl )-SAVE 
END IF 
65 CONTINUE 

63 CONTINUE 
YMAXU- YUK 1 6 0 ) -»-l 0 . 0 

XLOW--300 
XHI--295 
DO 80 J-1,120 
COUNT- 0 

DO 70 I-1,NSA«PLES 

IF (Y(I) .GE. XLOW .AND. Y( I ) .LT. XHI) THEN 
COUNT-COUNT+1 
END IF 
70 CONTINUE 

YV( J)-( FLOAT( COUNT )/NSAMPLES) *100 . 

XLOW-XLOW+5 
XHI-XHI-^5 
80 CONTINUE 

Calculating a scaling parameter for the vertical velocity histogram 

DO 82 1-1,120 
YVK I )-YV( I ) 

82 CONTINUE 
DO 83 1-1,120 

DO 85 J-1,119 

IF (YVKJ^l) .LT. YVKJ)) THEN 
SAVE- YVK J) 

YVK J )-YVl( J-»-l ) 

YVK J4l)-SAVE 
END If 
85 CONTINUE 

83 CONTINUE 
YHAXV-YVK 120) + 10. 0 

CALL HIST{YU,YV,YMAXU,YMAXV) 

C 

35 CLOSE(UNIT-ll ) 
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WRITE( * , * ) ' DO YOU DESIRE A PLOT OF THE COMPUTED VELOCITY', 

1 ' VALUES? (Y/N - DEF-Y)' 

READ( * , ' ( A1 ) ' ) ANS 

IF (ANS .EQ. 'N' .OR. ANS . EQ . 'n') THEN 

GO TO 30 
END IF 

CALL PT_PLOT( NSAMPLES, INFILE) 

WRITE( •, • ) 'DO YOU DESIRE TO VARY THE PLOTTED VELOCITY RANGE?', 
1 ' (Y/N — DEF-Y)' 

READ( *, ' ( A1 ) ' ) ANS 

IF (ANS .EQ. 'N' .OR. ANS .EQ. 'n') GO TO 30 
CALL PPLOT(NSAMPLES, INFILE) 

C • Again? 

30 WRITE!*,*)' ' 

WRITE (LUNTI,*) 'Again? (O-NO, 1-YES) ' 

READ (LUNTI,*) lYESNO 
IF ( lYESNO.EQ.l ) GO TO 100 
C 

RETURN 

C * Foraats; 

700 F0RMAT(5X,F9.4,1X,F9.4,8X,I4) 

800 F0RMAT(5X,F9.4,1X,F9.4) 

900 FORMAT!//, 2X, 'XPOS-' ,F6.2,2X, ' ZPOS- ' , F6 . 2 ) 

950 FORMAT! 2X, 'TEST SECTION VELOCITY- ', F5 . 1 , 2X , I 5 , ' SAMPLES') 

975 FORMAT! 2X, 'BLUE AND GREEN MIXING FREQUENCI ES— ' , 2 F5 . 1 , ' ( MHz ) ' ) 
END 



I 
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SUBROUTINE PT PLOT( NSArtPLES , INFI LE ) 



Written by S. D. Hedrick in July, 1987 to plot computed velocity values 
corresponding to airfoil phase angles. This subroutine is intended to 
provide the user with graphical information on acquired LDV velocities 
and assist in making a determination regarding the adequacy of the data 
obtained at the current probe position. 



IMPLICIT NONE 
INTEGER NSAMPLES,! 

REAL*4 XVAL( 4000 ) ,YVALU( 4000 ) ,YVALV( 4000) ,XCONTR( 4 ) ,YCONTR{ 4 ) 

CHARACTER*10 INfILE 

character*!! XLABEL 

CHARACTER* 15 YLABEL 

CHARACTER* 19 TITLE 

character*! dummy 

OPEN { UNIT-11 , FILE- INFILE, STATUS- 'OLD' ) 

DO 10 I-l, NSAMPLES 

READdl ,100) YVALU( I) , YVALV( I ) ,XVAL( I ) 

XVAL { I ) - XVAL ( I ) /1 0 . 0 
10 CONTINUE 

100 FORMAT(5X,F9. 4,1X,F9. 4,8X,F4.0) 

Set up the coordinate axis 

XLABEL-' PHASE ANGLE' 

YLABEL-' VELOCITY (ft/s)' 

TITLE-' COMPUTED VELOCITIES' 

XCONTR( 1 )-6 .0 
XCONTR{2)-0.0 
XCONTR(3)-360.0 
XCONTR( 4 ) -45.0 
YCONTRd )-6.0 
YCONTR(2)--200.0 
YCONTR{ 3)-600.0 
YCONTRi 4 )-100 .0 

CALL LGPSPLOTd, 'EXSY' , XVAL, YVALU,0, XLABEL, YLABEL, , , ,XCONTR, 

1 YCONTR, , TITLE) 

Plot the individual data points on the axis 

CALL LGPSPOiNTd, XVAL, YVALU, .03, 3, NSAMPLES, , , ) 

CALL LGPSPOINTd, XVAL, YVALV, .03, 4, NSAMPLES, , , ) 

CALL LGP$PUT_TEXT(1,4.5,5.75, '* U-COMPONENT' , , ) 

CALL LGP$PUT_TEXT( 1,4. 5,5.5, '0 V-COMPONENT' , , ) 

WRITE(*,*) 'Press <CR> to terminate plot' 
read(*,200) dummy 

200 rORMAT(Al) 

CALL LGPSTERMINATE PLOT(l,l) 

CLOSE(UNIT-ll) 

RETURN 

END 
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SUBROUTINE PPLOT( NSANPLES , INF I LE ) 

IMPLICIT NONE 
INTEGER NSAMPLES, I 
REAL VMIN,VMAX,INTV 

REAL*4 XVAL( 4 000 ) , YVALU( 4 000 ) , YVALV( 4 000 ) ,XCONTR{ 4 ) ,YCONTR( 4 ) 

CHARACTER*10 INFILE 

CHARACTER*!! XLABEL 

CHARACTER*! 5 YLABEL 

CHARACTER*!9 TITLE 

CHARACTER*! DUMMY 

OPEN {UNIT-1! , FILE- INFILE, STATUS- 'OLD' ) 

DO !0 I-!,NSAMPLES 

READ( !! ,!00) YVALU( I ) , YVALV( I ) ,XVAL( I ) 

XVAL( I )-XVAL( I )/!0.0 
10 CONTINUE 

100 FORMAT(5X,F9.4, !X,F9.4,8X, F4.0) 

Set up the coordinate axis 

WRITE(*,*) 'ENTER VELOCITY RANGE DESIRED (»in,Bax)' 

READ(*,*) VMIN,VMAX 
INTV-{ VMAX-VMIN)/5. 

XLABEL- 'PHASE ANGLE' 

YLABEL- 'VELOCITY (ft/s)' 

TITLE- 'COMPUTED VELOCITIES' 

XCONTR( 1 )-6.0 
XCONTR(2)-0.0 
XCONTRC 3 )-360.0 
XCONTR( 4 )-45.0 
YCONTR( ! )-6 . 0 
YCONTRC 2 )-VMIN 
YCONTR( 3 )-VMAX 
YCONTR( 4 )-INTV 

CALL LGP$PLOT{ !, 'EXSY' ,XVAL,YVALU, 0, XLABEL, YLABEL, , , ,XCONTR, 
! YCONTR, , TITLE) 

Plot the individual data points on the axis 

CALL LGP$POINT(!,XVAL,YVALU, .03,3, NSAMPLES, , , ) 

CALL LGP$POINT(!,XVAL,YVALV, .03, 4, NSAMPLES, , , ) 

CALL LGP$PUT_TEXT(1, .05, .3, 'ASTERISK-U COMPONENT', , ) 

CALL LGP$PUT_TEXT(1, .05, .1, 'CIRCLE-VXOMPONENT' , , ) 

WRITE(*,*) 'Press <CR> to terminate plot' 

READ(*,200) DUMMY 
200 FORMAT(Al) 

CALL LGP$TERMINATE_PLOT( 1 , 1 ) 

CLOSE(UNIT-l!) 

RETURN 

END 
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SUBROUTINE HIST( YU, YV, YMAXU , YM^XV) 



Written by S. D. Hedrick in August, 1987 to plot a histogram of the per- 
centage of data points which fall within specified velocity limits. The 
subroutine produces histograms for each of two velocity components. Scal- 
ing parameters for the vertical axis and data arrays are computed within 
a preceding subroutine (HWWA). 



IMPLICIT NONE 

REAL *4 YU(160),YV(120) ,XC0NTRU( 4 ) ,XCONTRV( 4 ) ,YCONTRU( 4 ) 

REAL* 4 YMAXU, YMAXV, YCONTRVI 4 ) 

REAL XLOW,XHI ,XLOWU( 160),XHIU(160) ,XLOWV( 120),XHIV(120) 

INTEGER I 

CHARACTER*29 UTITLE 
CHARACTER*27 VTITLE 
CHAilACTER*10 XLABEL 
CHARACTER*18 YLABEL 
CHARACTER*! DUMMY 

Setting up the coordinate axis 

XLABEL- 'VELOCITIES' 

YLABEL-' PERCENTAGE OF DATA' 

UTITLE-'HORIZONTAL VELOCITY HISTOGRAM' 

VTITLE-'VERTICAL VELOCITY HISTOGRAM' 

XCONTRU( 1 )-6 . 0 
XCONTRU( 2 )--200 . 0 
XCONTRU( 3 )-600 . 0 
XCONTRU( 4 )-100.0 
XCONTRV( 1 )-6 .0 
XCONTRV( 2 )--300 . 0 
XCONTRV( 3 )-300. 0 
XCONTRV( 4 )-100. 0 
YCONTRU( 1 )-6 .0 
YCONTRU( 2) -0.0 
YCONTRU( 3 )-YMAXU 
YCONTRU( 4 ) -10.0 
YCONTRV( 1 )-6 .0 
YCONTRV(2)-0.0 
YCONTRV( 3 )-YMAXV 
YCONTRV( 4 )-10.0 

CALL LGP$PLOT( 1 , ' IXSY' , , YU , 0 , XLABEL , YLABEL , , , ,XCONTRU, 

1 YCONTRU, , UTITLE) 

Setting up arrays for horizontal and vertical histogram bar parameters 

XLOW— 200.0 
XHI — 195.0 
DO 10 1-1,160 
XLOWU( I ) -XLOW 
XBIU( I )-XHI 
XLOW-XLOW+5.0 
XHI-XHI<fS.O 
10 CONTINUE 

XLOW--300.0 
XHI— 295.0 
DO 20 1-1,120 
XLOWV( I )-XLOW 
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XHIV( I )-XHI 
XLOW-XLOW-i-5 . 0 
XHI-XHI^5 . 0 
20 CONTINUE 

Plotting the horizontal velocity histogram 

CALL LGP$HIST( 1 ,XLOWU,XHIU, YU,160, , , 2 ) 

WRITE(*,*) 'PRESS <CR> TO TERMINATE PLOT' 

READ( * , ' ( Al) ' ) DUMMY 

CALL LGP$TERMINATE_PLOTU , 1 ) 

Repeating the procedures to plot the vertical velocity histogram 

CALL LGP$PLOT(l, 'IXSY' , , YV, 0 , XLABEL , YLABEL , , , ,XCONTRV, 

1 YCONTRV, ,VTITLE) 

CALL LGP$HIST( 1,XL0WV,XHIV, YV,120, , ,2) 

WRITE!*,*) 'PRESS <CR> TO TERMINATE PLOT' 

READ! * , ' ! Al ) ' ) DUMMY 

CALL LGP$TERMINATE_PLOT(l,l) 

RETURN 

END 
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PROGRAM ORGIZE 



Written by S. D. Hedrick in July, 1987 to sequentially organize data files 
to obtain a least value of phase angle (third column of data), and associ- 
ated velocities, at top of file. This rearrangement of data was desired 
to assist with subsequent manipulation and analysis. In September, 1967 
this routine was incorporated as the central program in the processing and 
analysis of unsteady LDV velocity data generated by an oscillating airfoil. 
The program calls assisting subroutines to disgard bad data points and pro- 
vide wake profile plots. ^ 



IMPLICIT NONE 

REAL X(4 000), Y( 4000) , XPOS , ZPOS , VELO, UVL , WL , UNDI M( 30) ,VNDIM( 30) 
REAL ZPO( 30) ,SAV1 ,SAV2,MENU,MENV, PREJ,MEANU( 18 ) ,MEANV( 18 ) 

REAL MPHAVG( 18) 

INTEGER Z( 4000) , NSAMPLES , PHASA , I , J , K , N , M , SAVE , RE J , STDY 
CHARACTER*10 INPILE 
CHARACTER*! ANS 
STDY-0 

WRITE( * , * ) ' IS ANALYSIS fOR NON-OSCILLATING AIRPOIL? ( Y/N ) ' 

READ( * , ' ( A1 ) ' ) ANS 

IF (ANS .EQ. 'Y' .OR. ANS . EQ . 'y') STDY-1 

WRITE( *,*) 'ENTER NUMBER OF DATA ACQUISITION POINTS IN PROFILE', 
1 ' ( 30 maximum ) ' 

READ(*,*) N 

WRITE( *,*)' ENTER PHASE ANGLE DESIRED FOR ANALYSIS (degrees)' 
READ(*,*) PHASA 
DO 70 K-1,N 

WRITE( *,*)' ENTER NAME OF INPUT DATA FILE (LDVXXX.DAT)' 

READ( *, ' ( AlO) ' ) INFILE 

WRITE(*,*) 'ENTER NUMBER OF SAMPLES IN DATA FILE' 

READ(*,*) NSAMPLES 

OPEN (UN IT- 11 , FILE- INFILE, STATUS- 'OLD' ) 

DO 10 I-l, NSAMPLES 

READ(11,100) X( I ) , Y( I ) ,Z( I ) 

10 CONTINUE 

READ(11,200) XPOS, ZPOS 
READ(11,300) VELO 
IF (STDY .EQ. 1) THEN 
CLOSE(UNIT-ll ) 

CALL STEADY(NSAMPLES, INFILE, MENU, MENV) 

UNDIM( K )-MENU/VELO 
VNDIM( K)-MENV/VELO 
GO TO 80 
END IF 

M-NSAMPLES-1 
DO 20 J-1, NSAMPLES 
DO 30 I-1,M 

IF (Z( I>1 ) .LT. Z( I ) ) THEN 
SAVE-Z( I ) 

SAVl-X(I) 

SAV2-Y(I) 

Z(I)-Z(I>1) 

X( I)-X(H.l) 

Y(I)-Y(I>1) 

2( !♦! )-SAVE 
X( !♦! )-SAVl 
Y( I4l )-SAV2 
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END IF 
30 CONTINUE 

20 CONTINUE 

CLOSEIUNIT-ll ) 

OPENIUNIT-15, FILE- INFILE, STATUS- 'OLD' ) 

DO 40 I-1,NSAMPLES 

WRITE(15,100) X(I),Y(I),Z(I) 

40 CONTINUE 

CLOSE(UNIT-15) 

WRITE( * , • ) 'DO YOU DESIRE A VELOCITY/PHASE PLOT OF 
1 ' DATA? (Y/N — DEF-Y)' 

READ( * , ' ( A1 ) ' ) ANS 

IF (ANS .EQ. 'N' .OR. ANS .EQ. 'n') GO TO 50 
CALL PT_PLOT(NSAMPLES, INFILE) 

50 CALL BDPTS( NSAMPLES, INFILE, MENU, MENV, RE J, PREJ,HEA; 

1 MPHAVG) 

CALL AVG(NSAMPLES, INFILE) 

WRITE( * , * ) 'DO YOU DESIRE INTERPOLATION OF THE DAT; 
READ( • , ' ( A1 ) ' ) ANS 

IF (ANS .EQ. 'N' .OR. ANS .EQ. 'n') GO TO 55 
CALL INTERP(NSAMPLES, INFILE) 

55 CALL CVFIT(NSAMPLES, INFILE, PHASA,UVL,WL, MENU, MEN\ 

1 MEANU,MEANV, MPHAVG) 

UNDIM( K )-UVL/VELO 
VNDIM( K)-WL/VELO 
80 ZPO(K)-2POS 

70 CONTINUE 

CALL PROPLOT ( UNDIM , VNDIM , ZPO , XPOS , N ) 

Format statements 

100 TORMAT( 5X, F9 . 4 , 1 X , F9 . 4 , 8X , I 4 ) 

200 FORMAT (//, 7X, F6 . 2 , 7X , F6 . 2 ) 

300 FORMAT( 24X,F5.1 ) 

STOP 

END 
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SUBROUTINE BDPTS ( NSArtPLES , INF I LE , MENU , MENV , REJ , PREJ , 

1 n EANU , M EANV , « P H A VG ) 

REAL X(4000),Y(4000), HEANU( 18 ) , MEANV( 18 ) , STDVU( 18 ) 

REAL STDW( 18 ) ,«PHAVG( 18 ) , MINU , MINV , ftAXU , MAXV , MENU , MENV 

INTEGER Z( 4000) , PHASEA , PHASEB , SUM , COUNT , NSAMPLES , NSAMP , REJ 

REAL* 4 XCONTR( 4 ) ,YCONTR( 4 ) ,PREJ 

CHARACTER*! ANS , DUMMY , DUMM I 

CHARACTER*10 INFILE 

CHARACTER*!! XLABEL 

CHARACTER*21 YLABEU 

CHARACTER*19 YLABEV 

CHARACTER*18 TITLE 

OPEN( UNIT-15, FILE- INFILE, STATUS- 'OLD' ) 

NSAMP-NSAMPLES 
DO 10 I-1,NSAMPLES 

READ(15,100) X(I),Y(I),Z(I) 

10 CONTINUE 

100 FORMAT ( 5X , F9 . 4 , IX , F9 . 4 , 8X , I 4 ) 

PHASEA-0 
PHASEB-200 
DO 20 1-1,18 
SUM-0 
SUMU-0. 

SUMV-0. 

PHAVG-0. 

DO 30 J-1,NSAMPLES 

IF (Z(J) .LE. PHASEB .AND. Z(J) .GE. PHASEA) THEN 
SUMU-SUMU^X( J) 

SUMV-SUMV+Y( J ) 

PHAVG-PHAVG+(FLOAT(Z( J) )/10. ) 

SUM-SUM4«1 

END IF 
30 CONTINUE 

IF (SUM .EO. 0) THEN 

MEANU( I )-0 . 

MEANV( I )-0. 

MPHAVG( I ) - ( FLOAT ( PHASEA)-t-100 . )/10. 

VARU-0. 

VARV-0 . 

GO TO 40 
END IF 

MEANU( I )-SUMU/FLOAT(SUM) 

MEANV( I )-SUMV/FLOAT( SUM) 

MPHAVG( I )-PHAVG/FLOAT( SUM) 

VARU-0 . 

VARV-0. 

DO 40 J-1,NSAMPLES 

IF (Z(J) .LE. PHASEB .AND. Z(J) .GE. PHASEA) THEN 
Al-X( J )-MEANU( I ) 

A2-Y( J)-MEANV( I ) 

Bl-Al**2/FLOAT( SUM) 

B2-A2**2/FLOAT( SUM) 

VARU-VARU^Bl 

VARV-VARV+B2 

END IF 
40 CONTINUE 

STDVU( I )-S0RT( VARU) 

STDW( I )-S0RT( VARV) 

PHAS EA- PHAS EA+ 200 
PHASEB-PHASEBo-ZOO 
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20 CONTINUE 

WRITE( * , * ) 'DO YOU DESIRE A STANDARD DEVIATION PLOT? ( Y/N - DEF-Y)' 
READ( *,200) ANS 
200 rORMAT(Al) 

IF (ANS ,EQ. 'N' .OR. ANS . EQ . 'n') THEN 

GO TO 50 
END ir 

Setting up the coordinate axis for the plot 

XLABEL-' PHASE ANGLE' 

YLABEU-' HORIZONTAL VELOCITIES' 

YLABEV-' VERTICAL VELOCITIES' 

TITLE-' STANDARD DEVIATION' 

XCONTR{ 1 )-6 . 0 
XCONTRC2 )-0.0 
XCONTRC 3 )-360.0 
XCONTR( 4 )-45.0 
YCONTR( 1 )-6 . 0 
YCONTRC 2 )--200. 0 
YCONTR(3)-600.0 
YCONTR( 4 ) -100. 0 

CALL LGP$PLOT( 1 , 'EXSY' ,MPHAVG,MEANU,0,XLABEL, YLABEU, , , ,XCONTR, 

1 YCONTR, , TITLE) 

Plotting the horizontal velocity data points on the axis 

CALL LGP$POINT(l,MPHAVG,MEANU, .03,4,18, , , ) 

Plotting the horizontal standard deviation 

CALL LGP$STNDEV{1,MPHAVG,MEANU,STDVU,18, ) 

45 WRITE(*,*) 'PRESS <CR> TO TERMINATE PLOT' 

READ(*,200) DUMMY 
CALL LGP$TERMINATE_PLOT{l,l ) 

Plotting the vertical standard deviation 

CALL LGP$PLOT(l, 'EXSY' ,MPHAVG,MEANV,0,XLABEL,YLABEV, , , ,XCONTR, 

1 YCONTR, , TITLE) 

CALL LGP$POINT( 1 ,MPHAVG,MEANV, . 03 , 4 , 18, , , ) 

CALL LGP$STNDEV( 1 , MPHAVG . MEANV , STDW , 1 8 , ) 

47 WRITE(*,*) 'PRESS <CR> TO TERMINATE PLOT' 

READ(*,200) DUMMI 
CALL LGP$TERMINATE_PLOT{ 1 , 1 ) 

C 

50 COUNT- 1 
PHASEA-0 
PHASEB-200 
DO 60 1-1,18 

MINU-MEANU(I)-2.0*STDVU{ I ) 

MINV-MEANV{ I )-2 . 0*STDW{ I ) 

MAXU-MEANU(I)+2.0*STDVU{ I ) 

MAXV-MEANV{I) + 2.0*STDW(I ) 

DO 70 J-1,NSAMPLES 

If (2{J) .LE. PHASEB .AND. Z(J) .GE. PHASEA) THEN 
If (X(J) .LT. MINU .OR. X(J) .GT. MAXU) THEN 
COUNT-COUNT^X 
IF ( J .EQ. 1 ) GO TO 70 
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DO 80 K-0,J-2 
X( J-K)-X( J-K-1) 

Y( J-K)-Y( J-K-1) 

2U-K)-2( J-K-1 ) 

80 CONTINUE 

GO TO 70 
END IF 

IF (Y(J) .LT. MINV .OR. Y(J) .GT. MAXV ) THEN 
COUNT-COUNT-^1 
IF (J .EQ. 1) GO TO 70 
DO 85 K-0,J-2 

J-K)-X< J-K-1) 

Y( J-K)-Y{ J-K-1) 

2{ J-K)-2( J-K-1) 

85 CONTINUE 

END If 
END IF 
70 CONTINUE 

PHASEA-PHASEA-f 200 
PHASEB-PHASEB+200 
60 CONTINUE 

CLOSED unit-15 ) 

OPEN ( UNIT-12 , FILE-INFILE, STATUS- 'OLD' ) 

DO 90 I-COUNT,NSA«PLES 

WRITE(12,100) X( I ) , Y( I ) ,Z( I ) 

90 CONTINUE 

NSANPLES-NSAMPLES+1 -COUNT 
CLOSED UNIT-12 ) 

OPEN < UNIT-1 3 , FILE-INFILE, STATUS- 'OLD' ) 

SUU-0. 

SUV-0. 

DO 95 I-1,NSAMPLES 

READ(13,100) X(I),Y(I),2(I) 

SUU-SUU-fX{ I ) 

SUV-SUV-^Y{ I ) 

95 CONTINUE 

MENU-SUU/FLOAT( NSAMPLES ) 

MENV-SUV/FLOAT( NSAMPLES ) 

CL0SE(UNIT-13 ) 

REJ-NSANP-NSAJ1PLES 

PREJ-( 1 FLOAT(NSAMPLES ) /FLOAT ( NSAMP) ) ) *100 . 

PHASEA-0 
PHASEB-200 
DO 96 1-1,18 
SUH-0 
SUMU-0. 

SUMV-0. 

PHAVG-0. 

DO 97 J-1 ,NSA?<PLES 

IF (2(J) .LE. PHASEB .AND. Z(J) .GE. PHASEA) THEN 
SUHU-SUMU-^X( J) 

SUMV-SUMV-^Y( J ) 

PHAVG-PHAVG+(FLOAT(Z( J) )/10. ) 

SUM-SUM+1 
END IF 

97 CONTINUE 

IF <SUH .EQ. 0) THEN 
HEANU( 1 )-0. 

NEANV( I )-0. 

«PHAVG( I )-(FLOAT(PHASEA)^100. )/10. 

» GO TO 98 

END IF 

MEANU( I )-SUMU/FLOAT( SUM) 

MEANV( I )-SUMV/FLOAT( SUM) 

MPHAVG( I )-PHAVG/FLOAT< SUM) 

98 PHASEA-PHASEA+200 
PHASEB-PHASEB+200 

96 CONTINUE 
RETURN 
END 
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SUBROUTINE AVG ( NSAMPLES , I N F I LE ) 

REAL X( 4000 ) , Y( 4000 ) 

INTEGER Z( 4000) , COUNT , SUM , NSAKPLES 
CHARACTER*10 INFILE 
COUNT- 1 
M-NSAMPLES-1 

OPEN( UNIT-12, FILE- INFILE, STATUS- 'OLD' ) 

DO 10 I-1,NSAMPLES 

READ(12,100) X( I ) ,Y( I ) , Z( I ) 

10 CONTINUE 

100 FORMAT( 5X,F9.4 ,1X,F9.4,8X, 14 ) 

DO 20 I-1,M 

IF (2(1+1) .EQ. Z(I) .AND. Z(I+2) . NE . Z(I)) THEN 

COUNT-COUNT+1 
X(I+l)-(X(I+l)+X(I))/2.0 
Y( I+l )-( Y( I+l )+Y( I ) )/2.0 
DO 30 J-0,I-2 
X( I-J )-X( I-J-1 ) 

Y(I-J)-Yd-J-l) 

Z( I-J )-Z( I-J-1 ) 

30 CONTINUE 

END If 

IF (Z(I+2) .EQ. 2(1)) THEN 
SUM-1 

DO 25 J-2,NSAMPLES-I 

IF ( Z( I+J) .EQ. Z( I ) ) THEN 
SUM-SUM+1 
END IF 
25 CONTINUE 

DO 27 J-1,SUM 

X( I + SUM)-X( I + SUM)+X( I + SUM-J) 

Y( I+SUM)-Y( I+SUM)+Y( I+SUM-J) 

27 CONTINUE 

X( I+SUM)-X( I+SUM)/( SUM+1 ) 

Y( I+SUM)-Y( I+SUM)/(SUM+1) 

DO 28 J-1,I-1 

X( I+SUM-J)-X( I-J) 

Y( I+SUM-J)-Y( I-J) 

2( I+SUM-J)-2( I-J) 

28 CONTINUE 
COUNT-COUNT+SUM 

END IF 
20 CONTINUE 

CLOSE( UNIT-12 ) 

OPEN(UNIT-14 , FILE- INFILE, STATUS- 'OLD' ) 

DO 40 I-COUNT,NSAMPLES 

WRITE(14,100) X( I ) , Y( I ) ,2( I ) 

40 CONTINUE 

NSAMPLES-NSAMPLES+1-COUNT 

CLOSE(UNIT-14) 

RETURN 

END 



) 
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SUBROUTINE I NTERP ( NSArtPLES , INFILE) 

IMPLICIT NONE 
CHARACTER*10 INFILE 
REAL X( 3600) ,Y( 3600) 

INTEGER 2(3600) , NSAMPLES , I , J , COUNT , NSA , FAC 
CHARACTER*! ANS 

OPEN( UNIT-U , FILE-INFILE, STATUS- 'OLD' ) 

DO 10 I-1,NSAMPLES 

READ(11,100) X(I),y(I),2(I) 

10 CONTINUE 

100 FORMAT( 5X,F9.4,1X,F9.4,8X,I4) 

COUNT-0 
NSA-NSAMPLES 
DO 20 1-1,3600 

IF (COUNT .GE. NSA) THEN 
2( I )-I-l 

X( I ).( (X( 1 )-X( I-l ) )/(2( 1 )-2( I-l ) ) ) *(2( I)-2(I-1 ) )+X( I-l ) 

Y( I ).( ( Y( 1 )-Y( I-l ) )/(2(l )-2( I-l ) ) )*( 2( I)-2( I-l ) )+Y( I-l ) 
NSAMPLES-NSAMPLES+1 
GO TO 20 
END IF 

IF (2(1) .EQ. I-l) COUNT-COUNT+1 
IF (2(1) .NE. I-l) THEN 
NSAMPLES-NSAMPLES+1 
DO 30 J-0,NSAMPLES-I-1 

X(NSAMPLES-J )-X(NSAMPLES-J-l ) 

Y( NSAHPLES-J )-Y(NSAMPLES-J-l ) 

2(NSAMPLES-J )-2(NSAMPLES-J-l ) 

30 ' CONTINUE 

2( I )-I-l 

IF (I .EQ. 1) THEN 

X(I)-((X(I-i-l )-X( NSAMPLES) ) /( 2 ( I ♦ 1 ) -2 ( NSAMPLES ) ) )• ( 2( I )- 
1 2 (NSAMPLES) ) +X ( NSAMPLES ) 

Y( I )-( ( Y( !*►! )-Y( NSAMPLES) )/( 2 ( I + l ) -2 ( NSAMPLES ) ) )*( 2( I )- 
1 2(NSAMPLES) ) +Y( NSAMPLES ) 

GO TO 20 
END IF 

X(I)-((X(l4l)-X(I-l))/(2(I+l)-2(I-l)))*(2(I)-2(I-l))+X(I-l) 
Y( I )-( ( Y( I + l )-Y( I-l ) )/(2( I + l )-2( I-l ) ) )* (2( I )-2( I-l ) )"^Y( I-l ) 
END IF 
20 CONTINUE 

CLOSE( UNIT-11 ) 

WRITE( *, * ) 'DO YOU DESIRE TO REDUCE THE NUMBER OF INTERPOLATED', 
1 ' DATA POINTS BY AN INTEGER FACTOR OF 3600? ( Y/N ) ' 

READ( * , ' (A1 ) ' ) ANS 

OPEN (UNIT-12, FI LE-INFILE, STATUS- 'OLD' ) 

IF (ANS .EQ. 'Y' .OR. ANS . EQ . 'y') THEN 
WRITE( *,*)' ENTER INTEGER FACTOR' 

READ(*,*) FAC 
FAC-20 

DO 34 I-FAC, NSAMPLES, FAC 

WRITE( 12, 100 ) X( I ) , Y( I ) ,2( I ) 

34 CONTINUE 

NSAMPLES- 3600/FAC 
GO TO 40 
END IF 

DO 40 I-I, NSAMPLES 

WRITE(12,100) X(I),Y(I),2(I) 

40 CONTINUE 

CLOSB(UNIT-12) 

» RETURN 

END 
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SUBROUTINE FOURI ER ( F , M , A , B , N ) 

DIMENSION F(0:M),FA(0:300),FB (0:300) 

DATA PI/3.14159265/ 

H-2»PI/M 
TYPE*, H 
J-0 

DO 2 I-0,M 
C-4 

IF (I .EQ. J) C-2 

IF (I .EQ. J) J-J*2 

FA( I )-F( I ) *COS(N*H*I ) 

FB( I )-F( I ) *SIN(N*H*I ) 

IF (I .EQ. M) T3-FA(I) 

IF (I .EQ. M) TT3-FB(I) 

IF (I .EQ. M-1) T2-FA(I) 

IF (I .EQ. M-1) TT2-FB(I) 

IF (I .EQ. M-2) Tl-FA(I) 

IF (I .EQ. M-2) TTl-FB(I) 

IF (1 .EQ. M-3) TO-FA(I) 

IF (1 .EQ. M-3) TTO-FB(I) 

IF (1 . EQ . 0 ) GO TO 2 
FA( I )-FA( I-l )*C*FA( 1 ) 

FB( I )-FB( I-l )*C*FB( I ) 

2 CONTINUE 
A-H/3./PI*(FA(M)-T3) 

B«H/3./PI*(FBlM)-TT3) 

IF { J .NE. M+1 ) GO TO 3 

A-H/3 ./PI*( FA(M)-4 . *T3-2 . *T2-4 . *T1-T0) 
A-A»3 . *H/8 ./PI*(T0*3 . *Tl+3. *T2+T3 1 
B-H/3 ./PI*( FB(M)-4 . *TT3-2 . *TT2-4 . *TT1-TT0) 
B-B+3 . *H/8./PI*( TTO+3 . *TT1*3. *TT2+TT3) 

3 CONTINUE 
RETURN 
END 
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SUBROUTINE CVFIT( NSArtPLES , INFILE, PHASA , UVL , WL , MENU , MENV , REJ , PREJ , 
1 MEANU , MEANV , MPHAVG ) 



Written by S. D. Hedrick in October, 1987 to fit a sinusoidal curve to LDV 
velocity data generated by an oscillating airfoil. The subroutine, utiliz- 
ing FOURIER, adapted from June 1977 OCS Airfoil program, fits a single 
cycle sine wave to the data and computes the u- and v-velocity components 
at a phase angle specified for the analysis. 



IMPLICIT NONE 
CHARACTER*! ANS , DUMMY 
CHARACTER*10 INFILE 

INTEGER NSAMPLES, I , PHASA , NSA , REJ , ARG 

REAL FU (0:3600), FV(0:3600),AU,AV,BU,BV, PI, X(3600),Y(3600) 

REAL UVL, WL, MENU, MENV, PREJ, MEANU ( 18 ) , MEANV ( 18 ) , MPHAVG ( 18 ) 

DATA PI/3.1415927/ 

WRITE( *,*)' DESIRE A VELOCITY/PHASE PLOT OF THE REFINED DATA?' 
WRITE ( *,*)'( Y/N--DEF-Y) ' 

READ( * , ' ( A1 ) ' ) ANS 

IF (ANS .EQ. 'N' .OR. ANS .EQ. 'n') GO TO 7 
CALL RP LOT (NSAMPLES, I N F I LE , RE J , PRE J ) 

7 OPEN{ UNIT-11 , FILE-INFILE, STATUS- 'OLD' ) 

DO 5 I-l, NSAMPLES 

READ( 11,100) X( I ) , Y( I ) 

5 CONTINUE 

100 FORMAT ( 5X, F9 . 4 , IX, F9 . 4 ) 

CLOSE(UNIT-ll ) 

DO 10 I-l , NSAMPLES 
FU( I-l )-X( I ) 

FV( I-l )-Y( I ) 

10 CONTINUE 

NSA-NSAMPLES-1 

CALL FOURIER( FU,NSA,AU,BU, 1 ) 

CALL FOURIER(FV,NSA,AV,BV,l) 

DO 20 I-0,NSA 

FU( I + l )-MENU-^AU*COS( ( 2 . * P I/FLOAT ( NSAMPLES ) ) ‘FLOAT ( I ) ) + 

1 BU*SIN( ( 2 . *Pl/FLOAT( NSAMPLES ) ) *FLOAT( I ) ) 

FV( I4l )-MENV-i«AV*COS( ( 2 . *PI/FLOAT( NSAMPLES ) ) ‘FLOAT ( I ) )♦ 

1 BV‘SIN( ( 2. ‘PI/FLOAT( NSAMPLES) )‘ FLOAT ( I ) ) 

20 CONTINUE 

WRITE( ‘ , ‘ )' PLOTS WILL NOW BE PROVIDED TO COMPARE SINUSOID WITH', 
1 ' VELOCITY DATA. PRESS <CR> TO CONTINUE.' 

READ( ‘ , ' ( A1 ) ' ) DUMMY 

CALL PLOT ( NSAMPLES, INFILE , fU , FV , MENU , MENV , MEANU , MEANV, 

1 MPHAVG) 

ARG-NSAMPLES* ( FLOAT( PHASA )/360. ) 

WRITE( ‘, ‘ ) ' ARG -' ,ARG 
UVL-FU( ARG) 

WL-FV( ARG) 

RETURN 

END 
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SUBROUTINE PLOT ( NSArtPLES , I NFI LE , FU , TV , MENU , MENV , MEANU , 
1 HEANV,MPHAVG) 



Written by S. D. Hedrick in October, 1987 to plot a comparison of velocity 
data and a sinusoidal curve fit. 



IMPLICIT NONE 

INTEGER NSAMPLES, I , BOUND , STOPU , STOPV 

REAL *4 XVAL( 4000 ) ,YVALU( 4000 ) ,YVALV( 4000 ) ,XCONTR( 4 ) , YCONTU( 4 ) 
REAL*4 YCONTV(4) , FU( NSAMPLES ) , FV( NSAMPLES ) 

REAL * 4 YMI NU , YMAXU , YM I NV , YMAXV , MENU , MENV 

REAL MEANU( 18 ) ,MEANV( 18 ) ,MPHAVG( 18 ) 

CHARACTER^IO INFILE 
CHARACTER*!! XLABEL 
CHARACTER*15 YLABEL 
CHARACTER*21 UTITLE 
CHARACTER*19 VTITLE 
CHARACTER*! DUMMY 

OPEN( UNIT-11 , FILE- INFILE, STATUS- 'OLD' ) 

DO 10 I-l, NSAMPLES 

READ( 11 , 100 ) YVALU( I ) , YVALV( I ) ,XVAL( I ) 

XVAL( I )-XVAL( I )/10 . 

10 CONTINUE 

100 FORMAT( 5X, F9 . 4 , IX, F9 . 4 ,8X, F4 .0 ) 

XLABEL- 'PHASE ANGLE' 

YLABEL- 'VELOCITY (ft/S)' 

UTITLE-'HORIZONTAL VELOCITIES' 

VTITLE- 'VERTICAL VELOCITIES' 

YMINU-MENU-150. 

YMAXU-MENU*^150 . 

YMINV-MENV-150 . 

YMAXV-MENV+150. 

XCONTR(1)-6.0 
XCONTR( 2 )-0 . 0 
XCONTR( 3 )-360.0 
XCONTR( 41-45.0 
YCONTU( D-6.0 
YCONTU( 2)-YMINU 
YCONTU( 3 ) -YMAXU 
YCONTU( 4 )-50.0 
YCONTV( D-6.0 
YCONTV( 2 )-YMINV 
YCONTV( 3 )- YMAXV 
YCONTV( 4 )-50.0 

CALL LGP$PLOT( 1, 'EXSY' ,XVAL,FU, NSAMPLES, XLABEL, YLABEL, , , , 

1 XCONTR,YCONTU, , UTITLE) 

CALL LGP$POINT(1,MPHAVG,MEANU,.03,3,18, , , ) 

WRITE( *,*)' PRESS <CR> TO TERMINATE PLOT' 

read(*,200) dummy 

CALL LGP$TERMINATE_PLOT( 1,1 ) 

CALL LGP$PLOT(l, 'EXSY' ,XVAL,FV, NSAMPLES, XLABEL, YLABEL, , , , 

1 XCONTR, YCONTV, , VTITLE) 

CALL LGP$POINT( 1,MPHAVG,MEANV, .03, 3,18, , , ) 

WRITE(*,*) 'Press <CR> to terminate plot' 
read(*,200) dummy 

200 FORMAT(Al) 

CALL LGPSTERMINATE PL0T(1,1) 

CLOSE<UNIT-ll ) 

RETURN 

END 
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SUBROUTINE PROPLOT ( UNDI M . VNDI Z PO , XPOS , N ) 

IMPLICIT NONE 

REAL *4 UNDIM ( 30 ) , VNDIM( 30 ) , ZPO( 30 ) ,XCONTU( 4 ) ,XCONTV( 4 ) 

REAL* 4 YCONTR( 4 ) ,UNDI( 3000 ) , VNDI (3000),ZP(3000) 

REAL XPOS 
INTEGER N,NBR 

CHARACTER*29 XLABEU , XLABEV 
CHARACTER*32 YLABEL 
CHARACTER*27 UTITLE 
CHARACTER*25 VTITLE 
CHARACTER * 1 DUMMY , ANS 
CHARACTER*9 XPO 
CHARACTER*4 XXPO 

XLABEU-' VELOCITY RATIO (UAMINFIN))' 

XLABEV- 'VELOCITY RATIO ( VA^ ( I NF I N ) ) ' 

YLABEL- ' PROBE POSITION (IN. FROM CENTER)' 

UTITLE- 'HORIZONTAL VELOCITY PROFILE' 

VTITLE- 'VERTICAL VELOCITY PROFILE' 

XPOS-XPOS/3.0 

WRITE(UNIT-XXPO, FMT- ' (F4 . 1 ) ' ) XPOS 

XPO-'X/C -'//XXPO 

XCONTU( 1 )-6 . 0 

XCONTU(2)-0.0 

XCONTUC 3 )-l . 2 

XCONTU( 4 )-0. 2 

XCONTVd )-6.0 

XCONTV( 2 )--0 . 7 

XCONTV( 3 )-0 .7 

XCONTV( 4 )-0 . 2 - 

YCONTR( 1 )-6 . 0 

YCONTR( 2 )--2. 0 

YCONTR( 3 )-2 . 0 

YCONTR(4)-0.5 

NBR-10*(N-1 )+l 

CALL LGP$PLOT( 1 , 'EXSY' , UNDI M , ZPO , 0 , XLABEU , YLABEL , 

1 , , ,XCONTU, YCONTR, , UTITLE) 

CALL LGP$POINT( 1 ,UNDIM, ZPO, . 03 , 3 ,N, , , ) 

CALL LGPSPUT^TEXTd, .4,5.5,XPO, , ) 

WRITE( *,*)' PRESS <CR> to terminate PLOT' 

READ( *, ' (Al ) ' ) DUMMY 

CALL LGP$TERMINATE_PLOT(l d ) 

CALL LGPSPLOTd, 'EXSY' , VNDIM ,ZPO , 0 , XLABEV , YLABEL , 

1 , , ,XCONTV, YCONTR, , VTITLE) 

CALL LGPSPOINTd , VNDIM, ZPO, . 03 , 3 , ) 

CALL LGP$PUT_TEXT(1, .4,5.5,XPO, , ) 

WRITE( *,*)' PRESS <CR> TO TERMINATE PLOT' 

READ( *, ' (Al ) ' ) DUMMY 

CALL LGP$TERMINATE_PLOT( 1 d ) 

WRITE( * , * ) 'DO YOU DESIRE A SPLINE INTERPOLATION CURVE FIT OF', 
1 ' THE DATA POINTS? ( Y/N ) ' 

READ( * , ' ( Al ) ' ) ANS 

IF (ANS .EO. 'N' .OR. ANS .EQ. 'n') GO TO 10 
CALL LGP$PLOT(l, 'EXSY' , UNDIM , ZPO , 0 , XLABEU , YLABEL , 

1 , , ,XCONTU, YCONTR, , UTITLE) 

CALL LGPSPOINTd ,UNDIM,ZPO, .03, 3,N, , , ) 

CALL LGPSSPLINE QHC ( UNDIM , ZPO , N , 1 0 , UND I , ZP , ) 

CALL LGP$PL0TM(T, 'EXSY' ,UNDI ,ZP,NBR, ) 

CALL LGP$PUT_TEXT(1, .4,5.5,XPO, , ) 

WRITE(*,*) 'PRESS <CR> TO TERMINATE PLOT' 

READ(*, ' (Al) ' ) DUMMY 

CALL LGP$TERMINATE_PLOT(l , 1 ) 

CALL LGP$PLOT( 1 , ' EXSY' , VNDI M , ZPO , 0 , XLABEV , YLABEL , 

1 , , ,XCONTV, YCONTR, , VTITLE) 

CALL LGP$POINT(l,VNDIM,ZPO, .03, 3,N, , , ) 

CALL LGPSSPLINE QHC ( VNDIM , Z PO , N , 1 0 , VNDI , ZP , ) 

CALL LGP$PL0TM(T, ' EXSY' ,VNDI , ZP,NBR, ) 

CALL LGP$PUT_TEXT( 1 , . 4 , 5 . 5,XPO, , ) 

WRITE( *,*) 'PRESS <CR> TO TERMINATE PLOT' 

READ( * , ' ( Al ) ' ) DUMMY 

CALL LGP$TERMINATE_PLOT(l , 1 ) 

10 RETURN 
END 
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